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Abstract: The WW scattering into gravitino and gaugino is here investigated in the 
broken phase, by using both gauge and mass eigenstates. Differently from what is obtained 
for unbroken gauge symmetry, we find in the scattering amplitudes new structures, which 
can lead to violation of unitarity above a certain scale. This happens because, in the 
annihilation diagram, the longitudinal degrees of freedom in the propagator of the gauge 
bosons disappear from the amplitude, by virtue of the SUGRA vertex. We show that the 
longitudinal polarizations of the on-shell W become strongly interacting in the high energy 
limit, and that the inclusion of diagrams with off-shell scalars of the MSSM does not cancel 
the divergences. 
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1. Introduction 

The gravitino G is the gauge field of local super symmetry, namely of Supergravity [1]. In 
the contest of Cosmology, such a particle plays a relevant role in several scenarios. After the 
spontaneous breaking of local supersymmetry, through the so-called super Higgs mechanism 
the gravitino obtains a mass that is proportional to the breaking scale. Accordingly, 
niQ depends on the particular model considered, and in principle it may range from the eV 
scale up to the TcV scale and beyond [2]. In general, gauge mediation predicts the gravitino 
to be the lightest supersymmetric particle, or LSP [3]. Then, if R-parity is conserved, it can 
be a very attractive candidate for Dark Matter. In other theories, for instance in gravity 
mediation, the gravitino is unstable and it has a lifetime which is usually longer than 
100 sec. This means that it decays after the beginning of the Big Bang Nucleosynthesis 
(BBN), affecting the abundances of the primordial light elements and eventually spoiling 
the success of the BBN [4]. Such problems have been extensively studied in the literature, 
setting different bounds on the masses of both stable and unstable gravitinos [5] . 

Inflation may solve these problems, since it dilutes enormously the gravitino abundance 
and provides with a more natural range for mg, between O(lMeV) and 0(100 GeV) [6]. 
The former constraints can then be relaxed and the cosmologically relevant gravitinos were 
produced during reheating, right after the end of inflation, mostly through hard 2 — > 2 
scattering processes of particles in the primordial thermal bath [7]. In Ref.[8], the authors 
calculated the gravitino production rate in supersymmetric QCD at high temperature, to 
the leading order in the gauge coupling. Ten hard 2^2 scatterings with a gravitino in 
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the final state were considered. The total contribution, with appropriate modifications due 
to the finite temperature of the thermal bath [9, 10], provides the collision term for the 
Boltzmann equation, and therefore an estimate of the gravitino number density. The same 
approach was then applied to the case of the electroweak interaction in the high-energy 
limit. By considering massless W bosons, relevant contributions to the total gravitino num- 
ber density were obtained [11]. By taking into account such results, the BBN constraints 
from gravitino production have been recently updated in [12], and an analytical procedure 
that is alternative to the numerical method was proposed in [13]. 

In this paper, we study for the first time what happens when the gravitinos are pro- 
duced at a centre of mass energy that is comparable to the EW scale, by assuming a non 
vanishing mass of the W bosons. In contrast to previous investigations [8, 11], and in some 
analogy with [14], we find that if the gauge bosons are massive, the squared amplitude of 
the WW scattering contains new terms which violate the unitarity above a certain scale, 
Eq.(2.15). Such quantities do not factorize any mass splitting which would vanish in the 
SUSY limit, contrary to what is expected [15]. In fact, in the annihilation diagram the 
longitudinal degrees of freedom in the W boson propagator disappear from the amplitude 
by effect of the supergravity vertex. After a comparison with the massless case, we show 
that the longitudinal modes of the W bosons become strongly interacting at high energies 
and the divergences hold at any centre of mass energy, as reported in Eq.(2.36). 

We now recall that broken supergravity is the effective limit of a more fundamental 
theory, since it is valid only below the SUSY breaking scale. Our calculations show that 
unitarity is broken at the same order, e.g. at ~ C(10^^ GeV) for tuq « 0(1 MeV). The 
result of this paper is thus consistent with the entire energy range allowed by SUGRA, 
therefore our study is phenomenologically motivated. 

In fact, the WW scattering can be observed at the LHC as a secondary process, for 
instance in gluon fusion [16]. From a more cosmological viewpoint, the fact that the result 
holds at any energy would be interesting in scenarios with both high and low reheating 
temperature Tr. While the former actually constitutes a rather standard background, 
Tr ^ 0(10^ GeV) is favoured in recent investigations on baryogenesis and Dark Matter 
[17]. Nevertheless, it would be interesting to study our result from a more formal point 
of view, as it seems to be a general feature of broken supergravity. This provides with an 
interesting theoretical perspective. 

The present article is organized as follows. In Section 2 we calculate the squared 
scattering amplitude of the process 

W + W^ — >W^ + G, (1.1) 

by using the basis of gauge eigenstates, in the case of massive spin 3/2 gravitino. We 
find some anomalies leading to violation of unitarity Eq.(2.15), whose origin is discussed 
in Section 2.1, where an argument concerning the high energy limit provides Eq.(2.36). 
Finally, at Section 2.2 wc calculate the effective limit by considering a gravitino mass 
which is much smaller than the SUSY mass splitting. The gravitino then becomes the spin 
1/2 goldstino, and contrary to the case with spin 3/2, no terms proportional to l/m~ are 
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generated in the squared amplitude. As we briefly discuss, this follows naturally from the 
structure of the effective Lagrangian. 

In Section 3 we consider the possibility of including scalars in the WW scattering. 
The only viable option is to add off-shell particles, but in principle this is forbidden by 
supergravity, which in the basis of gauge eigenstates does not admit the coupling gravitino- 
scalar-gaugino. We then use the formalism of mass insertions to add new contributions 
to the squared amplitude, up to the second order in the perturbation theory developed 
around a mass splitting. The result is that the scalars (we considered the neutral Higgs 
doublet in the MSSM) not only do not remove the divergences, but they generate also new 
anomalies. 

In Section 4 the basis of mass eigenstates is taken into consideration. The process is 
now the following: 

W+ + W- + G, (1.2) 

where Xq is the i-th neutralino in the MSSM. The mass eigenstates include the contributions 
of the photon and of the Z^, which are exchanged in the annihilation channels, so this is a 
physical process that can be observed at the LHC. Nevertheless, even in this basis we find 
the same quantities appearing in Eq.(2.15), since the longitudinal modes of the Z° vanish 
from the amplitude and confirm the general statement that was formulated at Section 2.1. 
Moreover, the inclusion of the neutral Higgs scalars does not lead to the cancellation of the 
divergences, in analogy with the result of Section 3. 

All the calculations in this paper have been performed at the tree level (because 
SUGRA is a non-renormalizable theory) with the program FORM [18]. 

2. WW scattering in supergravity in the broken phase 

In this section, we study the weak process in the basis of SU(2)l gauge eigenstates, in 
analogy with the scattering of two gluons computed in [8], where the group SU(3)c was 
considered. In the broken phase, the gauge bosons are massive and, as we will see in the 
following, this will make the difference with respect to the case of QCD. 

The scattering of two gluons g"' and g'' into a gravitino G and a gaugino (gluino) g'^: 

g-ik)+g\k')^rip')+Gip), (2.1) 

where a, b, c are indices of the SU(3)c algebra, corresponds to four irreducible Feynman 
diagrams, as it is shown in figure 1. 

We regard this process as the massless limit of (1.1), which instead consists of two 
bosons W"- and in the initial state, producing a gravitino and a wino: 

Wik) + W\k') W^{p') + G{p), (2.2) 

where a, b, c are now indices of the SU(2)l algebra. As it is put into evidence by the 
Feynman rules which are present in the Appendix (figure 13), the Supergravity interactions 
are universal. This means that the form of the interaction vertex between gravitino and 
gluon, or gravitino and W is exactly the same. The Feynman diagrams for the process 
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Figure 1: Gluons g° and scattering into a gravitino G and a gluino g'^. 

G '"^'"/i/i/i,. ^^G ^^'"^Aa 

11" 





W W 



Figure 2: The four diagrams which contribute to W + W'' — > W + G. 

(1.1), which are reported in figure 2, are indeed topologically identical to those in figure 1 
and we expect similar amplitudes. The matrix element M can be written as the sum of 
the four subamplitudes corresponding to the above diagrams: 



M = Ms + Mt + Mu + M^, 



(2.3) 



where Mg, Mt, contain the exchange of a particle in the channels corresponding to the 
Mandelstam variables s, t and u [19]. By using the Feynman rules summarised in Appendix 
A, the analytic expressions of the four separated contributions result as follows^: 

M, = 1^ [r?°^(A; - k'f + r]^''{2k' + - r]'"'{2k + k')f^] x 



Mu 



g^abc 



s — m 

1 



w 



4Mp t - m'L 



g^abc 



1 



4Mp u-m'i 



{r^{pm',l'']r{^-f + m^)rv's,{p')} eS(A;)e^(A:') 



w 



g^abc 



(2.4) 



Here Mp = (SttGn)'^^'^ = 2.43 x 10^^ GeV is the reduced Planck mass and g and Cabc srG, 
respectively, the coupling constant and the structure constants of the group SU(2)l. 



^As we will discuss immediately, it can be shown that the amplitudes for the gluon scattering (2.1) turn 
out to be exactly the same, modulo constant factors. 
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Regarding the spinors, V'^(p) is the wave function of a gravitino with four-momentum 
p and heUcity state s, represented by a four-component Majorana spinor, with spin 3/2 
and Lorentz index v'^,{p') denotes an antiwino with four-momentum p', hehcity state 
s' and SU(2)l index c in the final state, represented by a spin 1/2 Majorana spinor (we 
remand to [20] for the expansion in plane waves and for a detailed discussion about the 
concept of " antiparticle" for Majorana spinors). In the following, '0^(p) and Vg,{p') will be 
replaced respectively by ip^ and v'^ for simplicity of notation. 

Finally, e^(fe) and e^(fe') are the polarization vectors of two W bosons with, respec- 
tively, four-momenta k and k', SU(2)l indices a and b and Lorentz indices a and (3. 

The polarization sum of an antiwino with momentum p', group indices c and n, helicity 
state s' and mass in our normalization is [20]: 

Y,vt'{p'mP') = i^'-m^)5'^\ (2.5) 

s' 

whereas for a gravitino with momentum p, helicity state s and mass m^, the general 
expression for the spin sum is [21]: 



(2.6) 

When the energies are much larger than the gravitino mass, the full gravitino projector 
(2.6) can be approximated by [8]: 

E ^»^t{p) « -fn,. + Ip^- (2.7) 

o in- 

s G 

In the following, we will anyway use the full spin sum Eq.(2.6) instead of (2.7), since we 

are considering low energy scales which may be comparable with the gravitino mass. 

By using the Ward identities, we now rewrite the total amplitude in a form that clearly 
reveals the difference between massless and massive W scattering in supergravity. Eqs.(2.4) 
provide with the following Feynman amplitude M for the process (1.1): 



M = MJ el{k)el(k') = f-f^ [ri'^f^(k - k'Y + r]P''{2k' + A;)" - r/"'^(2A; + k'f] 



X 

1 r., .„ n .. 1 



s — m 



w 

- ^^2^ {^^^[^'^ 7^]7n^ - / + ^iy)7"«^} - V^Mb", i^]i''vAe%{k)e^p{k'). (2.8) 

As it was already discussed, the universality of supergravity implies that the Feynman 
diagrams and the amplitudes for the scattering of massive W are similar to those for 
massless W (or for gluons). It can be easily proven that the tensor density Af"^^ is indeed 
exactly the same in both cases, modulo the factor in (s — m^) that does not modify 
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the overall result^. The Ward identities for the case considered, where the gauge bosons in 
the initial state have four-momenta ka and k'^, can be written as [22]: 



where Sab and T^b are Lorentz scalars, which we call Goldstone amplitudes, as they reflect 
the degrees of freedom which are carried by the Goldstone bosons. Equations (2.9) allow to 
write the squared amplitude, summed over the polarizations of the initial and final states, 
in an interesting form. A straightforward calculation gives: 

E l-^l' = E IK^aW4(^')P = E.p.n \Kb'? - ZsmniS'^Xb)- (2-10) 
spin spin 

The scalar amplitudes have the following expression: 



4A 

which allows to recast Eq.(2.10) as: 



Sab = ^[ -i— ^ = r„,, (2.11) 

4Mp \ s-m^ I 



EiA^P = EK^>T-Ei^»^i'- (2-12) 

spin spin spin 

Accordingly, the spin-averaged squared matrix element of the scattering 

W + W'' — yW'' + G, 
can now be obtained by using (2.8) and averaging (2.12) on the spins of the initial states: 



_2 Isabel 2 
spin -"^ 



(2.13) 

where mj can be either m-w or m^, and / is a function of dimension \mass^\. 

The result that was obtained in SUSY QCD [8] for the scattering of two gluons (2.1) 

is: 

spin spin ^ \ G / 



spin spin 

hence, by comparison of (2.13) with (2.14), we find the relevant difference in the term 



-t{s + t) 



(2.15) 



^This property is anyway not intuitive, and it constitutes the main point of our analysis. It is discussed 
into details in subsection 2.1. 
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The above quantity contributes to the differential and to the total cross section of the 
process as follows: 

It is then clear that (2.15) can violate the unitarity of the theory above a certain centre of 
mass energy, which depends on the gravitino mass. This is rather similar to what happens 
in the Standard Model, if the scalars are not included in the WW scattering [19]. However, 
we will see in Section 2.1 that the above effect is not related to missing scalars, but to the 
structure of supergravity interactions. Nevertheless, due to the reduced Planck mass Mp in 
the denominator, this energy is of the same order of the SUSY breaking scale, meaning that 
the phenomenology of the process (1.1) is not dramatically constrained. A detailed study 
of the mechanism generating the term (2.15) is anyway required, and possible cancellations 
of this unexpected effect should be investigated. This is what we are going to address in the 
following. We begin by considering the squared amplitude written in the form (2.12), so 
that it is possible to compare it to the case of SUSY QCD. By computing the corresponding 
quantity for the gluon scattering (2.1): 

^ = ^ |M„f p - Y.iT'^'s:,) - Y.{S^'t:,) = ^ |M„f p - 2 J] (2.17) 

spin spin spin spin spin spin 

the key clement can be identified in the different contribution of the scalar amplitudes. It 
can be shown that M^jf is formally the same of Eq.(2.8), that (2.11) is replaced by 

gab ^ 9sf^ ^MiiLh!!: j = ^"^ (2.18) 

where /"^"^ are the structure constants of SU(3)c, and that the masses of the W bosons in 
the kinematics do not give any contribution which would cancel (2.15). It is straightforward 
to check that (2.17) follows from the sum over the polarizations of a massless gauge boson 
with four-momentum k [22]: 

E = - {na. - + 5-^^T + L. (2.19) 

Similarly, Eq.(2.12) is derived from the projector of a W boson with momentum k, mass 
iTiw and helicity state r: 

Y.'f\k)e'y^\k) = - - ^) ^T + L. (2.20) 

The structure of (2.12) and of (2.17) makes it possible to separate the two distinct contribu- 
tions of the above projectors to the squared amplitude of the process. Jospin \ -^'^ab\'^ comes 
from r/ajy (i.e. from T) while the second term (or L), which represents the longitudinal 
degrees of freedom of the W and of the gluon^, gives Y^gpi^ I'S'abl^- 

^We remark that the gluons do not have any proper longitudinal dofs, since they are massless. Such 
degrees of freedom are the Faddeev-Popov ghosts. 
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Transverse (T) and Longitudinal (L) 


mw = 


myv / 


LL 


nt{s+t) 


4-1 i-, 1 4-\ 

3m1 

, Li 


TL 


r,t{s+t) 


t{s + t) 
'im?~ 


LT 


ntis+t) 

G 


t[s + t) 
3m~ 


TT 






Total contribution 





t{s + t) 
3m~ 



Table 1: The contributions of the spin sums (2.19) and (2.20) to Y^spin l-^P- 



Alternatively to writing the squared amplitude in the form (2.12), and in order to 
check the calculation, one may act simply by "brute force" and substitute the spin sums 
(2.20) directly into the square of (2.8). The divergent terms which are generated by the 
interferences of the transverse and longitudinal components of the spin sums are reported 
in Table 1, for both massless and massive W. They lead respectively to (2.13) and to (2.14). 

By adopting two different procedures, we have thus found a clear indication that (2.15) 
is related to the longitudinal degrees of freedom of the W, which are related to the mass. 
This effect will be discussed into details in the next subsection. 



2.1 About the longitudinal polarizations of the W bosons, and their behaviour 
in the high energy limit 

The result Eq.(2.15) is here considered. First we show what can be the origin of the 
problem, then by considering a supercurrent argument we prove that the divergence (2.15) 
appears not only in the electroweak phase, but also at high energies. 

The amplitude of the s-channel, with an off-shell W, contains the product: 

{k + k')a{k + k' 



—I 



m 



w 



Vau + - 1) 



{V5^(p)[^ + r,7l7'^<'(p')}- (2.21) 



Prom the structure of the W boson propagator in the gauge. 



< A^{q)A^{-q) >=< A^{k+k')Aa{-k-k') >= 



or in particular, in the unitarity gauge (^ — > oo): 



m 



w L 



riau + (e - 1) 



{k + k')^{k + k'). 



< A^{k + k')A„{-k-k') >-- 



—I 



s — m 



w 



Vau + 



{k + k\{k + k% 



m 



w 



(2.22) 



(2.23) 



we see that, for both massless and massive gauge bosons, the part of the propagator which 
takes into account the longitudinal degrees of freedom of the particle: 

(k + k')„{k + k')^' 



m 



w I 



2 

w 



(2.24) 
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is cancelled identically when contracted with the SUGRA vertex 



ViCAl,^) = [1^ + 7^ (2.25) 

by virtue of the commutator. This means that the result is the same in the Feynman 
(^ = 1), Landau (^ = 0) and unitarity gauge, and that in supergravity this happens for 
every diagram with an off-shell gauge boson, either massless or massive. This fact has no 
consequences in the former case, since the longitudinal degrees of freedom are unphysical 
and cancel anyway during the calculation^. On the contrary, if the gauge bosons are 
massive such cancellations should not occur, as the longitudinal polarizations of the W 
are physical. It is then expectable that in our case this mechanism would generate in the 
squared amplitude terms such as 

spin G G 

which violates unitarity above a certain scale. As it is well known from the Standard 
Model, the Feynman and the unitarity gauge are not equivalent, unless we add diagrams 
with the exchange of a scalar to the amplitude written in the Feynman gauge [19]. This 
way, we recover the longitudinal modes that are missing if the propagator is expressed only 
as 

< AJk + k')Aa(-k - k') >= -i '^""^ , (2.27) 

and the divergences in the squared amplitude vanish, restoring unitarity. 

In contrast, the problem persists in SUGRA, since such a theory does not admit the 
coupling gravitino-scalar-gaugino in the basis of gauge eigenstates. As it will be shown 
in Section 3, to invoke the mixing of higgsinos and winos seems to be the only way to 
introduce the scalars in the process. Nevertheless, we will also show that the inclusion of 
the Higgs doublet of the MSSM does not cancel the anomalies. Clearly, we should recover 
the missing degrees of freedom and restore unitarity in some other way. 

Starting from general considerations, we can now discuss about the high energy be- 
haviour of the result (2.13). In SUSY theories with unbroken gauge symmetry, the super- 
current Sn is conserved in the SUSY limit. Namely, the right hand side of 

p^S^ = A(M), (2.28) 

where p'* is the four momentum of the gravitino and A(M) is the mass splitting of the 
supermultiplet, vanishes when all the masses are set to zero. However, in the electroweak 
theory the W mass cannot be set to vanish. In other words, we expect that in the process 
(1.1) the longitudinal modes of the W become strongly interacting in the limit niw 0, 
therefore in this case the statement (2.28) docs no longer hold. We will now prove this 
claim, which implies that our result holds at any centre of mass energy. 



''This it is well known from the Standard Model [19]. It was explicitly shown in [8] that the same holds 
also in supergravity. 
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The matrix element can be written as 
therefore the squared amphtude can be recast as 
When summed over the gravitino polarizations, it becomes 



(2.29) 



tp—spin 



ip—spin 



(2.30) 



Since a term with the gravitino mass squared in the denominator can appear only from the 
longitudinal spin 1/2 components of G in cither (2.6) or (2.7), one can replace the gravitino 
projector with the second term of Eq.(2.7), that is: 



"^.PuPfM 

oP 2^ 

3^ m~ 

G , 



(2.31) 



Therefore, it is always proportional to p^Sf^ = M.{ipn p^), i.e. to the matrix element 
with the replacement p^- 

It is now possible to demonstrate that the scalar density PuS^^ is not proportional to 
any SUSY breaking parameter, and accordingly that it does not vanish in the SUSY limit. 
In fact, by explicit calculation, one can show that Eqs.(2.4) imply the following: 



P^S^ 



g^abc 



1 



4Mp (s - m^){t - mi^){-s - t + 2m^ + m^) 
where the tensor density can explicitly be written as: 



(2.32) 



+st{-r]'^^$ - 7/°^^^' - 2A;^7° + 2A;'°7^ + 7°^7^|^ - 7°^7''|^' - 2p'°^7^ - 2p'^7°)] . 

(2.33) 

We have used equations like k"ea{k) = k'^efj{k') = 0, p'v'^{p') = m^v'^{p') = 0{m), and 
k'^ = k'^ = ruy^. At this stage, one could naively regard the terms proportional to as 
unimportant and make them vanish together with m^. The claim of the conservation of 
the supercurrent would immediately follow. 

However, this is not the case since, as we will see in a moment, the longitudinal 
polarization vectors e'f^{k) of the W bosons play a crucial role, as inferred in the discussion 
which follows Eq.(2.15). According to [19], we can write: 



mw 



niw 



(2.34) 
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which in the approximation A; — ^ oo can be recast, component by component, as 

e^W-^ + of^V (2.35) 



The components of are growing as k, in agreement with the general relations ■ k = 
and k ■ k = m^. It is then possible to replace the polarization vectors in (2.32) with the 
above expression^. By using the kinematics we then obtain, modulo constant factors: 



p^S^ _ (^ 



m 



v^ + 0[ , (2.36) 

where now m is really any mass. We can conclude that in the scattering of massive gauge 
bosons in SUGRA, the supercurrent is not conserved in the SUSY limit, and that the non 
vanishing term in the right hand side can contribute to generate the divergences which 
were found in the previous section. Eq.(2.36) completes the proof and constitutes the main 
result of this paper together with (2.15). 

We now open a parenthesis about the Ward identities, recalling what has been found 
in the previous section. As clearly described in Ref.[22], in a generic SU(2) gauge theory 
with fermions in a doublet representation, the unitarity of the iS-matrix implies certain 
conditions which any scattering amplitude has to satisfy. In particular, with the notations 
we have used in this section, the imaginary part of M^^ connects the initial and final states 
to all physical states with the same energy-momentum of the initial and final states. 

To summarize, in this section we have shown first that the structure of the supergravity 
vertex cancels the longitudinal degrees of freedom of any propagating gauge boson. Thus, 
if it is massive, terms with bad high energy behaviour should be generated. With an 
argument based on the supercurrent, we have then proven that the above result holds both 
at low scales and in the regime of high energies, where the longitudinal degrees of freedom 
of the W become strongly interacting. 

2.2 The effective Lagrangian and the limit of light gravitino 

After computing the scattering of massive W bosons into massive gravitinos, we now con- 
sider the effective limit of the process (1.1). In this approximation, the gravitino mass 
mQ is much smaller than the mass difference between bosons and fermions in the chiral 
and gauge multiplets, i.e. than the SUSY mass splitting in the observable sector [21]. By 
consequence, the helicity ±1/2 components of the gravitino field become dominant, and it 
can be regarded as a spin 1/2 fermion (the goldstino). The wave function is now rewritten 
as follows: 

V'^-n/^— (2.37) 
V 3 uIq 

where ip is the spin 1/2 goldstino. The spin sum changes from (2.6) to: 

Y^^^^){p)^i^)(p)=f±mQ, (2.38) 



^The corrections 0{m/E\^) provide with terms proportional to m^, which are absorbed into 
O {m? /mo). Namely, the proof given above has validity also at low energies. 
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where, as usual, the positive sign refers to fermions, while the negative sign refers to 
antifermions. The interaction Lagrangian is modified accordingly [8]: 

m~ — ml _ _ qm7^, 

(2.39) 

therefore the vertices are formally the same as in the Lagrangian of the full theory. It is 
easy to see that the above expression factorizes every goldstino vertex with a mass ratio 
which reminds the supersymmetric mass splitting. This results evident from the form of 
the spin-averaged squared amplitude: 



Spin ^ 




17^ 3 3t^\ ^ 
48' + 4^+47^+^^"'^^^^"'*^ 



(2.40) 



where mj is any mass and /(s, t) is a dimensionless function of s and t. The above result is 
consistent with (2.38) and with the approximation Eq.(2.39) of the Lagrangian. We have 
thus confirmed that, as expected in view of the goldstino spin sum (2.38), the divergences 
in the form (2.15) are generated only by the terms proportional to l/m~ in the projector 
(2.6), namely by the longitudinal spin 1/2 modes of the gravitino. 

3. Inclusion of scalars and mass insertions 

Up to now we have seen that, in the scattering of massive W bosons into massive gravitinos, 
one finds in the cross section (2.16) terms which violate unitarity, Eq.(2.15). It was shown 
in Section 2.1 that this happens because the longitudinal degrees of freedom of the W in 
the propagator vanish from the amplitude, by virtue of the supergravity vertex, and that 
the on-shell W becomes strongly interacting at high energies. 

Similarly to the Standard Model, the additional longitudinal modes of the W (which 
here are missing) may be taken into account by adding diagrams which contain scalar 
particles. In the case of (1.1), the only possibility is given by an annihilation channel 
with an off-shell scalar, but in the basis of gauge eigenstates there may be some formal 
difficulties. First of all, as it can be inferred from figure 2, the scalar (j)'^ has to be included 
in the adjoint representation of the gauge group, but in general (and in particular in 
Supergravity) the scalars belong to the fundamental representation. In this case, however. 
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the gauge symmetry is broken and the mixing of gauge eigenstates may allow the particle 
multiplets to pass from one representation to the other. We can then consider the dominant 
couplings of the gravitino with the supercurrent in the supergravity Lagrangian [1]: 



8Mp 



V5M[7^7^]7^A('^^i^i?, (3.1) 



where Xlr~ Pl,rX^ ~ 2 ^ anticipated in Section 2.1, the coupling of the 

gravitino to a scalar and a gaugino cannot be found in the above equation, but we can use 
the method of mass insertions and invoke the mixing of the gaugino with a chiral fermion 
that couples to the gravitino and to a scalar, as shown in figure 4. We now observe that 
SUGRA is a chiral theory, meaning that matter fermions are represented by spinors which 
are left-handed Weyl, not Majorana. This means that they have definite helicity, therefore 
the couplings in the interaction Lagrangian contain only one chirality projector, either 
or Pr, and provide with chiral vertices (figure 15 in the Appendix). On the other hand, 
the supergravity couplings in (2.4) regard Majorana fermions which do not have definite 
helicity. Namely, in the interference of the Higgs scalars with the original amplitude (2.8) 
that we now label with M^^\ there may occur cancellation of some polarizations which 
in principle would contribute to the process. Accordingly, we expect that the diagram in 
figure 4 would not cancel the divergences (2.15). 

Within the particle content of the MSSM, the 
role of the scalars in the WW scattering (1.1) can 
be played by the two neutral Higgs bosons and 
-fff, which break the gauge symmetry by taking the 
vacuum expectation values vi and V2- Let M^^^ 
be the amplitude that is the sum of two distinct 
diagrams with the exchange of a Higgs, as in figure 
3 (in figure 3 and in figure 4, the higgsino hi is 
labelled by a solid line, as we want to stress the 
fact that it is different from the wino). We must 
be also careful about the vertex of the Higgs with 
the W bosons, as it is not the same coupling of the 
Standard Model, since now we consider a doublet of scalars. This leads to the following 
expression for the vertex [23]: 





Figure 4: The diagram with the ex- 
change of a scalar 



1^1 



v,Z}l + V2Z-fl. 



2i 



(3.2) 



Zr is the rotation matrix transforming the gauge eigenstates and into the mass 
eigenstates H and /i, which are respectively the heavy and light Higgs in the MSSM: 




cos Q — sm Q 
sin a cos a 




(3.3) 
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We are actually interested only in the gauge eigenstates, therefore we impose ck = in 
(3.2), and obtain the following: 

. 2 

V{HfW^W^) = ^-e"-^%ir]'^^ , (3.4) 



P> 2 

The mixing factors can be found in the neutralino mass matrix [20]: 

Ai 

A2 



+ — tor hi, 



9V2 , r 
-— forhs. 



(3.5) 

(3.6) 
(3.7) 



where g is the coupling constant of the weak interaction. Moreover, the higgsino propa- 
gating in the external leg of each diagram provides with the factor 

1 1 



m 



(i = l,2). 



(3.8) 



w 



The above approximation will give a clearer form of the amplitude, and it can be safely 
assumed since the higgsino mass does not enter the kinematics and is not relevant for the 
final result. We thus obtain the amplitudes corresponding to the neutral Higgs bosons of 
the MSSM in the basis of gauge eigenstates: 



g^abc 






')7^1 - 




8V2Mp 




s 






g^abc 






i'h'^il - 




8V2Mp 











which by recalling the fundamental relation [20] 

4 

can be recast as a single amplitude: 



2 9 / 2 , 2\ 



(3.9) 
(3.10) 

(3.11) 



g^abc 

8V2Mp 



m 



w 



s + [g'^imjj^ - mjjjvl - 4m^m^J 



(3.12) 



It can be proven that adding the above expression to (2.8) does not remove the original 
anomaly in (2.15), since the dominant divergences result as: 



^2 Isabel 2 



3ml 



Ml 

s2(7s + 3t) 64s(s + 2t) Ss^ ' 

"5. 



64 X 48 \w?~ml m^~is + t) V2m^m^ 



(3.13) 
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W 




Figure 5: The t— and u— channel diagrams with the mass insertions. 



clearly showing that the interferences of M^^^ with M^^^ give only the term with \/2. 
However, by iterating the method of mass insertions, one may add further diagrams, devel- 
oping a perturbation theory around a mass splitting. At the first order, the matrix element 
is given by: 

m(i) =Mt,_, +M„,,,, (3.14) 

where Mt^ ^ ^^'^ 2 contain, respectively, the t and u channel. Each of them splits into 
two subdiagrams, since the neutral higgsinos hi and /12 both mix with the wino W ^ as it 
is shown in figure 5. The amplitudes result as follows: 



My 



*1,2 



1 g<^ahc ( 



1 



m 



+ 



t — m~ (t — m~)t 



4V2 Mp 

for the diagrams with the t-channel. The sum of the u-channels corresponds to: 



7V|6S(A.)6j(fe'), 
(3.15) 



M 



1 9(-ahc ( m 



"1,2 



4\/2 Mp \m 



''W 



"W 



r 



m 



+ 



u — m~ (u — m~)u 
w ^ w' 



^v^ eS(fe)6^(fc') 



(3.16) 

At the second order, the amplitude M^^^^ contains the double mass insertion with a gaugino 
in the external leg, corresponding to the diagrams which are shown in figure 6, where 
i = 1,2. The corresponding amplitudes are: 

1 / T r. W n „ „s 1 



s — m 



w 



(^^ [fi + 0, 7.] 7^^^^ + {V^M 7"] 7^1^' - / + m^h^v^} 



w 



u — 



{V^^[|^',7/3]7^(^-/ + m^)7-t;'=} -4'^^[7",7/^]7A'7;Me^(/c)e^(fcO, (3.17) 
'w / 

that is nothing but Eq.(2.8), multiplied by the factor m^^/m—. The total amplitude can 
be finally written as the sum of the four contributions: 

MiW" + — >W^ + G) = m(°) + m(^) + m(^) + m(2). (3.18) 

It can be shown that the squares and the interference of M^^^ and M^^^ give subdominant 
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Figure 6: The diagrams with two mass insertions in the external legs. 



terms which are proportional to e.g. 

f 'mw\ t{s + t) 



m~ I 3m~ 
^ w/ G 



(3.19) 



and that the leading divergences appear only in Eq.(3.13) and in 



y |M« + M(^)|2 ^ ^''^"'y 5'^^^^ • (3-20) 
' 48M2 V2m^m^ ^ ^ 

spin f 1 W G 

Interestingly, the above factor cancels exactly the corresponding term in Eq.(3.13). This 
means that the contribution of the amplitude M^^^ is suppressed, and that the Higgs 
diagrams contribute only with [M^^^p. We can then conclude that, as expected, the 

inclusion of the scalars does not restore unitarity since the cancellation of the longitudinal 
modes of the off-shell W boson still occurs, and the on-shell W bosons become strongly 
interacting at high energies. 

To summarize, up to this point we have considered a scattering process in Supergravity 
in the basis of gauge eigenstates of the weak interaction in the MSSM. It turns out that the 
squared amplitude contains terms which violate unitarity above a certain scale, Eq.(2.15). 
In the following section, we will see what happens if we consider instead the basis of 
electroweak mass eigenstates in the MSSM. 



4. WW scattering in the basis of mass eigenstates 

We now discuss an alternative approach with respect to the previous sections. In the 
Standard Model, for an unbroken gauge symmetry such as SU(3)c, the particles do not 
mix, but when the symmetry is (spontaneously) broken, the mass matrices become non 
diagonal. The mass eigenstates, new physical particles which can be detected at collider 
experiments, are then introduced when the mass matrices are diagonalized. This happens 
also in supersymmetric theories, and in particular in the MSSM. We will now discuss the 
WW scattering into gravitinos in such a basis, therefore the particle content now includes 
the four neutralinos x'i (^ runs from 1 to 4), the four cliarginos xf U = 1) 2) and the scalar 
Higgs multiplet of the MSSM. The process: 

W+ + W-^X°i+G, (4.1) 
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Figure 7: WW scattering in supergravity in tlie basis of mass eigenstates. 

where Xi is a neutralino with i running from 1 to 4, is the counterpart of Eq.(l.l) and 
corresponds to the set of diagrams in figure 7. From the above remarks, it is clear that this 
is a physical process that can be observed in experiments such as the LHC. First of all, 
we notice that the supergravity vertices in Appendix A now factorize the mixing factors of 
the mass eigenstates, but they do not change analytically. Namely, we expect to confirm 
the results obtained at Sections 2 and 3, due to the cancellation of the longitudinal modes 
in the Z boson propagator. More into details, with the notations of [2] and [20], we choose 
the following basis of neutral mass eigenstates: 

iPf = i-iX.,,-iX,,XH,Xj,), (4.2) 

that contains the photino 7, the zino z and the supersymmetric partners of the heavy 
H and light h Higgs scalars, the higgsinos H and h . The four particles mix, and the 
non-diagonal mass matrix Y' can be diagonalized by a unitary matrix A^': 

X° = Nl^i^f, N'*Y'N'-^ = Nd, (4.3) 

with i,j running from 1 to 4. In the following, the mixing factors will be labelled respec- 
tively by A^,-, A^;-, A^'- and A^^ , where the index i labels the ith neutralino. The charginos 

7' 2'' Hi hi 

can be instead defined as: 

xt = V,,^|^|, x7 = Ui,i^;, (4.4) 

where i = 1,2 and U, V are unitary matrices which are chosen so that 

U*XV-^ = Md, (4.5) 

with Md being a 2x2 diagonal matrix, with non negative entries. The chosen basis contains 
four fermions, which are made up of gaugino and higgsino components and mix in analogy 
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Figure 8: The diagrams which represent the amphtudes Ms and M^- 



with the neutral mass eigenstates, namely 

i:^ = {W+,H+,W-,Hj). (4.6) 

In the following, we will label with V^^y+ and U_^^_ the mixing elements of, respectively, 
the ith positive chargino with the wino and of the ith negative chargino with the wino 
W~ . We first consider the exchange of a gauge boson and the four-particle vertex, where Aj 
is the wave function of the i-th neutralino (figure 8). By taking into account Eq.(4.3) and 
the Feynman rules reported in the Appendix, we can write the corresponding amplitudes 
as follows: 



4:Mp 



gcosgty ^^a^^j^ _ ^ ^/3-(2A;' + kT - r/"'^(2A; + k'f] x 



AMp 



1 



rrir 



9 



(4.7) 

(4.8) 
(4.9) 



It is then clear that in the longitudinal modes in the Z propagator are cancelled by 
the SUGRA vertex, with exactly the same mechanism that was discussed in the previous 
sections. 

The novelty in Eq.(4.9) is the term N:^., i.e. the mixing factor of the neutral wino W"^ 
with the zth neutralino, in the basis given by 

{B,W^hi,h2), (4.10) 

which is alternative to (4.2), with the bino B. Both and B are linear combinations of 
7 and z, and the unitary matrix N diagonalizes the mass matrix Y, like in Eq.(4.3). 
My and can be recast into a single amplitude Mg. 



[rj'^^ik - k'Y + r/^"(2A:' + k)^ - <"(2fc + k'f] {\a^ + J^\^„] i:^] x 



cos ow 



m^N'~^ sin 9w 



s(cos2 Ows — w? 



w' 



(4.11) 
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Figure 9: The diagrams corresponding to the amphtudes Mt and M„ 



where mz = mw / cos Ow has been also used. As far as the t and u channels in figure 9 are 
regarded, the amplitudes result as follows: 



a 



v.. 



8Mp j^^t 



mi. I 



Of:-(l-75)+Of^(l + 75) 



9 



U. 



8Mp 



u — m 



mi L 



o,^*(i + 75) + o,; 



75) 



(4.12) 



(4.13) 



The heavy H and light h CP-even Higgs bosons are included into the diagrams in figure 
10, and correspond to the amplitudes: 



Mh- 
Mh 



2y/2Mp 
gmw 



N'-.cos(p-a) 

Hi ^ ' 

"tV^^ Sin (/?-«)" 



2 ^ 



m 



H 



2V2Mp 

which can be recast into the following: 
gmw 



A,(l-75)7n^ + r)V'. ^ 

1-^ ' 



mt 



(4.14) 
(4.15) 



Mh 



2y/2Mp 



N'^. sin (/? - a)(s - m|^) + N'^. cos (/? - q)(s - m|) 

A.(l - 75)7^1^ + r)V^. ./3, , 
(s-m^)(s-m^) 



(4.16) 



The Feynman amplitude of the process + — > + G can then be written as: 

M=Ms+Mt + M^ + M^+ Mmggs = M^tux + Mmggs- (4.17) 

We are now interested in checking whether divergences such as (2.15) appear also in this 
case. Accordingly, in the following we will not list the finite terms, but consider only the 
dominant anomalies. It is evident that the spin-summed square of the above object is very 
complicated, as it depends on N^^, Ofj and similar quantities. As the mixing factors have 
not yet been determined experimentally, we are anyway free to fix some constraints on these 
parameters. First we report the squared amplitude of the four diagrams corresponding to 
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+ 




Figure 10: The diagrams with exchange of the Higgs bosons. 

Mstux, then we will consider also the contribution of the Higgs bosons. By using the 
following equation: 



N^. = Nl^i sin Ow + N'^i cos Bw-, 



'71 



(4.18) 



we find cancellation of all the divergences which are proportional to l/m^m~ and to 
I/ui'Iy. Moreover, in the case which corresponds to 



Wi' 



2N^^{0^* + Of) + (lOi^f + lO^f ) = 0, 



(4.19) 



also the terms which are proportional to l/rn^f and to l/m^m~ vanish. If instead U. 
and are taken to be different, Eq.(4.19) is rewritten as: 



w- 



u 



jw- 



' 1 



(4.20) 



and no anomalies are cancelled. In any case, the terms which factorize l/m~ still remain, 
and the result is the following: 

Esp^n\Ms + Mt + M^ + M.,\^ 



9 



Ml cos2 Ow 3m| 



7« 



cos^ 9w 



-2cos2 6w 



2\N^.\His + t) + 2iV~ (Of + Of ).2 - 40,§Of } , (4.21 



where we considered the case of Eq.(4.19) for simplicity. It is straightforward to check that 
Eq.(2.15) is obtained in the limit when 9w — > and N:^.,Ofj,Ofj e°-'"^, namely in the 
basis of SU(2)l gauge eigenstates. We remark that during this calculation no mass has 
been set to zero, and that no identities other than Eqs.(4.18) and (4.19) have been adopted. 

The diagrams with the propagating Higgs bosons, corresponding to figure 10 and to 
the amplitude (4.16), contribute to the leading divergences only with the square: 



Nlsm{/3-a)+N'^.cos i/3 - a)\^ x 



Ml 12^^' hi 

2s^ 
cos^ Owm?- 



+ 



(2cos2 Ow - 1) 



(4.22) 
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as the interferences of Mniggs with the other four amphtudes provide only subleading 
contributions. The following equation [20]: 



Am^m^ cos^ 



2(5 



(4.23) 



together with mz = mw/ cosOw, has been also used. We have then found some analogy 
with Eqs.(3.13) and (3.20), namely the scalars contribute to the factor and reintro- 

duce the divergences proportional to l/rjiyy and to l/m^m~. It can be explicitly shown 
that the general case Ujy^^_ 7^ ^jw+ '^^^^ change the result (this is evident anyway, 
since (4.22) depends only on s). Eqs.(4.21) and (4.22) complete the calculation in the basis 
of the mass eigenstates, and they confirm what has been obtained in Section 2. 

One can also check whether the Ward identities (2.9), and accordingly gauge invariance, 
hold in the present situation as well. In the basis of mass eigenstates, suitable constraints 
on the mixing factors eliminate the unwanted terms in the contractions of and k'^^ with 
Mgtux in (4-17). The scalar density 5' is thus found in analogy with the results of Section 2. 
However, the Higgs bosons complicate things once again, as it is clear from the following: 



J k^W^ = Sk'^ + Snk^ 
I k'pW^ = Sk"" + Snk"^ ' 

where the scalar amplitudes S and Sh correspond to: 



(4.24) 



S = 



Sh 



gmw 
V2Mp 



AMp s — m^^ 
N'-.sm. {P - a)(s - m%) + N'-. cos (/3 - a){s - m^) 



(4.25) 

A,(1-75)V'mP'^ 
{s-ml){s-m%)' 
(4.26) 



In Eq.(4.24) there is a clear violation of gauge invariance and therefore of unitarity, by 
effect of Sh- Interestingly, the contraction of the longitudinal modes of the in the 
propagator, that is 



—I 



m 



w I 



{k + k')^{k + k% 



with the 3-bosons vertex 

'rfl^{k - k'Y + 77^^ (2A;' + k^ - rf^ilk + k'f 



(4.27) 



(4.28) 



if non vanishing, would have provided with additional terms proportional to k'^^ and k^ (or 
to A;" and /c'"), thus contributing to both the scalar densities S and Sh- The fact that the 
Higgs amplitudes cause violation of gauge invariance and of unitarity, confirms our result 
(4.22) and is consistent with the discussion in Section 3. 

To summarize, we have found that the squared amplitude of the scattering 



w+ + w- 
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contains terms which violate the unitarity of the theory above a certain scale. This happens 
because the longitudinal polarizations of the TP in the propagator vanish from the amplitude 
(4.8), and the Higgs bosons do not cancel the divergences. As it is discussed in Section 2.1, 
the contribution of the longitudinal polarizations of the on-shell W bosons is determinant, 
since in the high energy limit they become strongly interacting. In this section we have 
then shown that it is not possible to restore unitarity by introducing the suitable scalar 
particles which are found in the spectrum of the MSSM. 

5. Conclusions and outlook 

In this paper, we have studied in detail the WW scattering in Supergravity in the broken 
phase, and compared our results with the case of massless gauge bosons. If in the high- 
energy limit the W bosons are considered massless a priori, as in Ref. [11], the scattering 
amplitude does not contain any problematic terms, resulting in a formal analogy with the 
case of QCD [8]. 

In Section 2 it is shown that, in contrast, at low energies and with massive W, we find 
new structures, which could lead to the violation of unitarity in both the bases of gauge 
and mass eigenstates. This happens because in the annihilation diagram, the longitudinal 
degrees of freedom of the W boson in the propagator disappear from the amplitude by 
virtue of the supergravity vertex, implying that the longitudinal polarizations of the on- 
shell W become strongly interacting at high energies. This is discussed in Section 2.1, 
where we derive equation (2.36) that constitutes the main result of this paper, together 
with (2.15). In Section 3 it is then shown that, in the basis of gauge eigenstates, no scalar 
particle can provide with terms which would cancel (2.15). The same happens in the basis 
of mass eigenstates of the weak interaction in the MSSM, as discussed in Section 4. It is 
found indeed that the doublet of neutral Higgs bosons does not restore unitarity. 

As remarked in the Introduction, the result of this paper is however phenomenologically 
interesting, since SUGRA is the effective limit of a more fundamental theory and is valid 
only below the SUSY breaking scale, that here coincides with the energy at which unitarity 
is violated. Accordingly, we can apply Eqs.(2.13), (4.21) and (4.22), together with the other 
finite contributions of the WW scattering, to collider Physics and to Cosmology. In fact, 
the process we have analyzed in this article can be observed at the LHC as a secondary 
reaction, for instance through gluon fusion. Moreover, from a cosmological viewpoint, our 
result contributes to the thermal gravitino production in the early universe both at low and 
high reheating temperatures, as we have proven in Section 2.1 that it holds at any scale. 
Accordingly, it may be interesting to consider what happens also in other processes with 
on-shell W bosons, and to calculate their contribution as well. All these considerations 
constitute the background for future research. 

Nevertheless, from a more formal point of view, wc would expect to find that the 
divergence (2.15) is cancelled by some mechanism, which docs not seem to be as immediate 
as the inclusion of a scalar. This theoretical perspective provides the motivation for future 
investigations. 
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Appendix 

A. Feynman rules for the MSSM and supergravity 

The relevant Feynman rules used in our calculations are here summarised. The propagators 
in the gauge and the vertices of Supergravity and of the Minimal Supersymmetric 
Standard Model follow Refs. [21], [19] and [20]. 

In the following, Mp = (SvrGAr)-^/^ = 2.43 x 10^^ GeV is the reduced Planck mass, g 
and /"^^ are, respectively, the coupling constant and the structure constants of a generic 
gauge group. In our notations, Pl,r = ~ (l T 7^)- We point out that our expressions differ 
from those of Ref.[20] by an imaginary factor i, since the squared amplitude must be real. 
This is justified by the different normalization of the Lagrangians of Supergravity and of 
the MSSM. All the momenta are considered to run into the vertex and the fermion flow is 
denoted by a thin line with an arrow. 



Propagators and interaction vertices 



AAA/WWW 



5' 



ab 



k'^—m'i 



5^ 



3 I ^""^ 

Figure 11: The propagators of a vector boson with mass niA and of a scalar witli mass m^. 



a,M b,i^ k'^-m^j i j k'^-rrij^ 

Figure 12: The propagators of a Majorana fermion with mass tum and of a Dirac fermion with 
mass m£). 
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Figure 13: Gauge boson-gaugino-gravitino. 




r 4 r 

Figure 14: Gauge boson-gauge boson-gaugino-gravitino. 




■{PLYl''Vp)hj 



Figure 15: Scalar- Weyl fermion-gravitino. 





Figure 16: Gauge boson-gaugino-gaugino. 




-g^alO^PL + OffPnl 



Figure 17: Gauge boson-chargino-neutralino. 
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